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Abstract 

We perform one-zone simulations of the infall epoch of a pre-supernova stellar core in the pres- 
ence of neutrino flavor changing scattering interactions. Our calculations give a self-consistent 
assessment of the relationship between flavor changing rates and the reduction in electron fraction 
and re-distribution of initial electron lepton number among the neutrino flavors. We discuss and 
include in our calculations sub-nuclear density medium corrections for flavor changing scattering 
coherence factors. We hnd that flavor changing couplings e > 3 x 10^^ in either the <-> or 
i/g <-> channels result in a dynamically significant reduction in core electron fraction relatively 
soon after neutrino trapping and well before the core reaches nuclear matter density. 
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I. INTRODUCTION 



Core collapse supernovae are exquisitely sensitive to lepton number violating processes. 
This is because the infall (collapse) epoch of the pre-supernova core is characterized by low 
entropy l| and large lepton (electron and electron neutrino) degeneracy. Nearly all of the 
pressure support stems from these degenerate leptons. The effects of including neutrino 
flavor changing neutral current (FCNC) interactions in the infall stage of a core collapse 
supernova have recently been investigated in Ref. |2i]. It was noted there that neutrinos in 
the core of a collapsing star could undergo large numbers of scatterings due to the coherent 
amplification of the neutrino-quark flavor changing neutral current cross section for elastic 
scattering on heavy nuclei. Such interactions could cause significant numbers of electron 
neutrinos in the core to be converted to mu and tau neutrinos. In turn, this would open phase 
space for further electron capture and thereby significantly impact the pressure, homologous 
core mass, and the initial shock energy. 

The explosion of core collapse (Type II, lb, and Ic) supernovae is believed to be the result 
of gravitational collapse, subsequent hydro dynamic bounce of the star's core, and release of 
gravitational binding energy into neutrinos which ultimately provide the energy to revive 
and sustain the shock . One important feature of the model is that the 

entropy of the core is low [s/k ~ 1) and nucleons remain bound in nuclei during most of 
the collapse. The number of electrons in the core (hence, the pressure and homologous core 
mass) is governed by the electron capture reaction e~ + p ^ + n. When the neutrino 
mean free path becomes smaller than the size of the core (because of scattering on heavy 
nuclei) the neutrinos become trapped. They thermalize quickly and comprise a degenerate 
Fermi-Dirac sea. When the z/g Fermi level becomes high enough, electron capture is blocked 
and net reduction in (where Yj = {rif — nj)/ni,) no longer occurs on dynamical time 
scales. However, re-distribution of electron lepton number between i/g's and electrons will 
still occur as the density rises and the nuclear composition changes. 

Any further changes in the core's electron fraction during the collapse could result in a 
change in the collapse dynamics and explosion mechanism lOt] . Including neutrino FCNC 
interactions in the collapse model causes greater reduction in Ye during infall. This is because 
when electron neutrinos change fiavor by scattering, holes open in the z/g sea and the electron 
capture reaction can procede. 
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Neutrino-quark FCNCs of the form 

^ = ^^^7^^^g7.(4, + 4.7')g (1) 

were considered in Ref. Here, the parameters ey,^ and e\.. quantify the strength of 

the FCNC relative to the Fermi constant Gp- Current experimental constraints [n] on the 
FCNC couplings are ey^ < 10 for the channel Ue i^^ and ey^^ < 5 x 10 ^ for the 
channel Ue ^ Vt- (Similar, and in some cases better, constraints on these interactions may 
be possible from solar and atmospheric neutrinos [l^.) 

The cross section for neutrino flavor changing elastic scattering on heavy nuclei, mediated 
by the FCNCs of Eq. ([1]), was calculated in Ref. [2I and a coherent amplification was found. 
Using this cross section, and employingvalues of the coupling constant up to and beyond 
current experimental constraints, Ref. ^ gave estimates for the number of neutrino flavor 
changing scattering events which could occur in the core. The resulting reduction in Ye and 
implications for the stellar collapse model were then discussed in a qualitative sense. 

In this paper we present results of a one-zone calculation of the infall epoch of a pre- 
supernova star with neutrino-quark FCNCs included. Our code gives a more accurate ac- 

n 

counting of scattering rates and the change in Y^. than do the estimates of Ref. [j] and we are 
able to account for some of the feedback in the system. We model neutrino scattering with 
nuclei in the core medium and account for sub-nuclear matter density structure effects. By 
contrast, Ref. employed neutrino-nucleus vacuum cross sections with no accounting for 
medium effects. Reference {2! estimated which values of e would give a fast enough FCNC 
scattering rate such that reduction in Y^. would be possible. Here we actually compute what 
the reduction in Ye is for various values of e, including values below the best experimental 
bounds. We have discovered that maximal reduction in Ye is possible for values of e smaller 
than the best experimental bound in the Ve ^ channel, and that dynamically significant 
reduction in Ye is possible for values of e smaller than the best experimental bound in the 
^ channel. In section II we describe our code and method of computing the change 
in electron fraction. In section HI we discuss our results and their meaning for the stellar 
collapse model. In section IV we list the key approximations in our calculation and give 
an assessment of the possible impact of the potential uncertainties introduced by these. In 
section V we give conclusions. 
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II. ONE-ZONE CORE COLLAPSE SIMULATION 



We seek a self-consistent relationship between FCNC rates and the possible reduction 
in core electron fraction resulting from these processes. We simulate the core collapse with 
a one-zone calculation which computes reactions rates (including FCNC rates), thermody- 
namic quantities, equation of state (EOS) quantities, and electron and neutrino fractions. 
Though one-zone calculations obviously do not include a sophisticated treatment of hydro- 
dynamics or neutrino transport, and can contain many assumptions, they have been used 
successfully to model feedback between weak interactions and nuclear equation of state pa- 
rameters in the infall epoch of stellar collapse [isl . The validity they have rests on three 
key and non-controversial features of the infall epoch: low entropy; e~ and z/g degenerate 
conditions; and lepton capture rates dominated by the energetics scales associated with the 
high lepton Fermi levels. Note that these key features are also confirmed by large sophisti- 
cated numerical simulations 3, ^, Q] • 

A. Description of Calculation 

The code is a modified version of that used in Ref. \x3\. In the calculations done here 
and in Ref. (l3|, a single zone (with initial electron fraction Y^, density p, temperature T, 
entropy per baryon S, neutron mass fraction X„, neutron kinetic chemical potential fin, 
and neutron-proton kinetic chemical potential difference fi) is evolved assuming a uniform 
collapse rate. A standard Newton-Raphson algorithm is employed. As the density increases, 
the electron Fermi energy rises and the electron capture rate increases. At each density step 
the electron capture rate and collapse rate are used to find the change in electron fraction 
AYe, and then AY^ is used to estimate a change in entropy AS. The updated values for p, Y^, 
and S, along with explicit expressions for S and Xn are used by a routine which increments 
the temperature and utilizes a two-dimensional root finder to iteratively compute A„ and T. 
The mean nuclear mass A, fin, and p are also found during this iterative process. The EOS 
formulae used [isl for the mean nuclear mass, nucleon chemical potentials and nucleon- 
to-baryon ratios are based on a finite temperature liquid drop model (see Ref.s with 
a representative mean heavy nucleus and a sea of dripped neutrons. This is discussed in 
Appendix |Al 
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The core's electron fraction changes because of electron capture reactions. After neutrinos 
have become trapped in the core and the z/g's build up a degenerate Fermi sea, an equilibrium 
situation obtains: Yg and the net number of u^s per baryon Y^^ no longer change appreciably 
even though electron capture reactions, and the inverse reactions, are taking place. In the 
presence of FCNCs, u^s change flavor. As a result, phase space is opened allowing net 
electron capture to occur and causing further reduction in Yg and in overall electron lepton 
number. 

The physical reason for the reduction of Y^ is that electron capture reactions lower the 
number of electrons in the core. We count the reduction of Ye in two ways. The first way 
uses the electron capture rate for reactions occurring before equilibrium is established. The 
second way counts electron captures which occur as a result of phase space opening in the 
z/e sea, secondary to flavor changing scattering events. This is computed using the neutrino 
flavor changing rate. Counting the reduction in Yg in the first way applies only until complete 
beta-equilibrium obtains, while counting in the second way applies both before and after 
beta-equilbrium is established. We will discuss the first way here, and discuss the second 
way in Subsection III CI after we present the neutrino flavor changing scattering rate. 

The rates of electron capture on free protons and heavy nuclei are derived in Ref. [Ic 
and denoted respectively by Afp and Ah- The total rate of electron capture per baryon is 

^ = -X.Afp - ^A„, (2) 

where Xh ~ (1 — X„ — Xp) is the mass fraction of heavy nuclei and X„ and Xp are the 
neutron and proton mass fractions, respectively. (The number abundance of heavy nuclei 
relative to baryons is Yh = X^/ A while the corresponding abundances of the free nucleons 
are Yn = Xn and Yp = Xp.) We take the alpha particle mass fraction to be negligible, 
consistent with the low entropy infall conditions. Combining Eq. (1A12P for the collapse rate 
and Eq. ([2]) we have 



dYe 
dp 



1 - Xn - Xp 

— ApAfp — Ah 
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where pio = p/(10^°g/cm'^). This gives AY^ at each density step and is used to find Y^ until 
the beta equihbrium condition is imposed. 

At the onset of collapse, electron neutrinos created from electron capture stream freely out 
of the core. As nuclei become more neutron-rich and the cross section for ordinary coherent 



neutral current scattering becomes appreciable, high energy neutrinos begin to be trapped 
in the core and start to equilibrate. This occurs for a matter density of p ~ lO^^g/cm^. We 
start neutrino trapping at the density ptrap = 5 x lO^^g/cm^. When the simulation reaches 
density ptrap it begins accounting for the z/g's getting trapped in the core. The v^. fraction 
Yy^ is found by calculating AYg with equation ^ at each density step and imposing the 
condition AYJ^^ = AYg. This of course is an approximation because it means that after 
the "trapping density" is reached, every electron capture creates a neutrino which becomes 
trapped. In reality, neutrinos are not trapped instantaneously at some density, but rather 
are gradually trapped as neutrino diffusion times increase. Also, Vf. cross sections scale as 
neutrino energy squared so some lower energy neutrinos will still be escaping after the higher 
energy neutrinos have become trapped. Though our model for neutrino trapping is obviously 
siniplistic, it gives values for Yy^ consistent with the currently accepted core collapse model 

i3i,y,B,y,fly,3- 

One-zone collapse calculation results for various cases are shown in Tables Hl lIIII These 
give density pio, electron fraction "Kg, mu plus tau neutrino fraction Yy^ +Yy^, temperature T 
in MeV, entropy per baryon s/k in units of Boltzmann's constant, free neutron fraction Xn, 
neutron kinetic chemical potential p„ in MeV, mean nuclear mass A, mean nuclear radius 
Tnuc in fin, separation distance between nuclei Dsep in fm, average neutrino energy Ey in MeV, 
and average neutrino de Broglie wavelength Xy in fm. Note that our crude neutrino trapping 
and neutrino sea filling schemes give some unphysical results. For example, once neutrino 
trapping is enforced at ptrap = 5 x lO^^g/cm^, the entropy is taken as constant. However, 
the system is still out of chemical equilibrium so increments in density give a temporary 
and small drop in temperature. This gives a negligible overestimate of the FCNC effects 
near trapping because it produces slightly larger nuclei, with of order a few extra nucleons. 
Likewise, because we do not consider z/g-capture self consistently with electron capture, and 
because our liquid drop equation of state is inaccurate for high density and high neutron 
excess, runs with values e ^ 10^^ acquire a positive neutron kinetic chemical potential at 
sub-nuclear but large densities. Note however, that where this happens FCNCs have already 
had a significant effect. 

When the electron neutrino fraction reaches Yy^ = 0.05, we impose the condition that 
beta-equilibrium has been reached. We have chosen Yy^ = 0.05 as final equilibrium value of 
Yy^ to be consistent with the currently accepted supernova model and large scale numerical 
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simulations . At this point, in the absence of neutrino FCNCs and density- 

and composition-driven equihbrium shifts, the values for and Y^^ are the final values. 
This is another approximation since equilibrium does not actually obtain instantaneously. 
In Table [I] we have given results for a run of our simulation with flavor changing interactions 
turned off. The run used to produce this table started from a density of 3.7 x lO^g/cm^ 
and went to a final density of 3.8 x lO^^g/cm^. This table shows the changes in Y^ and Y^^. 
As soon as the simulation starts. Ye is decreasing. After the trapping density, Y^^ starts to 
increase and once it reaches 0.05, Yg and Y^^ no longer change. 



B. Neutrino Flavor Changing Rate in the Core 

Understanding the nuclear composition and equation of state in the core of a collapsing 
star is an active area of research. It is believed that as the core approaches nuclear matter 
density, p ~ lO^^g/cm^, the nuclear component undergoes a series of phase transitions as 
the individual nuclei merge and, in fact, eventually cease to exist [15]. During these phases, 
i.e. "pasta phases," the nuclear matter may take the form of rods, sheets or tubes. Recent 
work has focused on how neutrinos scatter coherently on stuctures in these phases. (See for 
example Ref. [3].) As outlined above, we use a liquid drop model to describe the nuclear 
component in the core. We use this model to describe the core only up to a density of 
p = 3.8 X lO^^g/cm^. The liquid drop model may not be valid over the whole density range 
where we have used it. By only running our simulation up to a maximum density which is 
an order of magnitude below nuclear density, we avoid most of the density range where it is 
guaranteed to be inaccurate. 

In Tables Hl lIIII we show values for the mean nuclear mass A, radius of the mean nucleus 
Tsep, and the nuclear separation distance, -Dgep- If should be kept in mind that a nuclear 
statistical equilibrium mix of nuclear sizes and masses will exist in the core. The mean 
nuclear mass is taken from Eq. flA4l) and the nuclear radius is found from Vnuc ~ A^^^ fm. 
To calculate the separation distance between nuclei, we assume each nucleus is in a Wigner- 
Seitz cell with cell volume Vc = l/riH, where uh = pN^Y^ is the number density of heavy 
nuclei and Na is Avagadro's number. Then Dgcp = "iRc, where Rc is the radius of the cell. 

In Table [H for example, we see that values for mean nuclear mass become as high as 
A ^ 300. At earlier epochs and lower densities in our simulation, we see that the nuclear 
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a foundation for modern full-scale supernova simulations [17]. At the higher densities, the 
nuclei have radii as large as 7 fm, and are separated by at least 37 fm. In the middle of 
the density range, the nuclei have radii less than 5 fm and are separated by at least 40 fm 
and as much as 80 fm. According to our liquid drop model, the values for A, rnuc, and D^ep 
indicate that we are not close to densities where the nuclei merge. Therefore, up to a density 
of 3.8 X lO^^g/cm^, it is reasonable to consider coherent scattering of neutrinos on individual 
nuclei in the core. However, in reality one should consider the more complicated problem of 
neutrino coherent scattering on the three dimensional structures in the pasta phases [igI ]. 

We calculate the rates for neutrino flavor changing scattering on free nucleons and on the 
mean nucleus. Coherent scattering on nuclei is the dominant flavor changing reaction. In 
the core it can be necessary to take account of scattering interference effects arising from 
conditions where more than one nucleus resides within a neutrino DeBroglie wavelength. 
This is an issue whenever 1/E,y ^ Rc [l8^, where E^, is the average neutrino energy. This 
condition means that the neutrino DeBroglie wavelength is comparable to or larger than 
the distance between nuclei. In a relativistically degenerate Fermi gas the average neutrino 
energy is 3/4 of the neutrino chemical potential /Xi,^, where ~ ll.lMeV{2piQY,yJ^^^ . In 
Tables Hl lIIII we show average neutrino energies and average neutrino DeBroglie wavelengths, 
X,y. Recalling that Rc = Dsep/2, and comparing Rc to X^, we see that we are in a regime 
where interference will occur. In Appendix |B] we present the cross section for neutrino 
scattering with nuclei in a medium. 

The neutrino flavor changing scattering rate employed here is as follows. First, the 
electron neutrino flavor changing scattering rate per mean nucleus is 



Here, p is the matter density, c is the speed of light and X is a factor which corrects for 
interference. This rate was obtained by multiplying the Uc flux by the cross section given in 
Eq. (1B5|) . The coherent ampliflcation factor is (2A^ + Z)^, where N and Z are the number 
of neutrons and protons in the mean nucleus. We have dropped sub and superscripts on the 
FCNC coupling e so that this rate is generic and refers to scattering on a rf-quark in either 
the Uc ^ or z/g ^ Vr channel. The rate for electron neutrino flavor changing scattering 




(4) 
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on free nucleons is 

\,^ = [pNaY,:)c'^2EI. (5) 

TT 

We do not concern ourselves with accounting for the fact that the free nucleons are in a 
medium; the cross section for coherent scattering on nuclei is larger by a few orders of 
magnitude, so the nuclei dominate the FCNC opacity. The total neutrino flavor changing 
rate per baryon is 



i 



A 2 



1 — Xn 



(6) 



C. Counting the Reduction in Electron Fraction 

Given that the parameter e is currently constrained by experiment to be smaller than 
~ 10~^ (z/g ^ z/t) or ~ 10^'^ iVf, ^ u^), the flavor changing scattering rates could be less 
than the electron capture rate. We work in a limit such that whenever an electron neutrino 
undergoes a flavor changing scattering and opens a hole in the z/g sea, the hole is immediately 
filled by a z/g produced via an electron capture. We argue later in this section that our 
calculated scattering rates justify this approximation. In this limit, FCNC transformation 
of electron neutrinos into mu and tau neutrinos will not change Y^^. Rather, the change in 
electron fraction AY^. is simply minus the total change in the sum of mu and tau neutrino 
fractions 

AFe = -A(n^ + nj. (7) 

In other words, the net reduction in Ye equals the sum of the net increase of the mu and tau 
neutrino fraction. The change in Y^^ + Y^^ in a density step is 

10-^2 / \ 

Pio Vg/cmV 

where pio = p/(10^°g/cm'^) and is given in Eq. ([6]). (This expression is analogous to 
Eq. O) The total reduction in Ye stemming from neutrino FCNCs is found from Eq.s 
and (IHl) by summing the increments from each density step. We impose the condition that 
the FCNC interactions are turned on in our one-zone simulation at the neutrino trapping 
density, which we take to be ptrap = 5 x lO^^g/cm^, a value consistent with large scale 
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simulations. Before trapping, neutrinos are freely streaming out of the core and net electron 
capture is not yet blocked. It would not matter to the core's final value of Ye if electron 
neutrinos changed flavor before streaming out of the core. 

In Figure [T] we have plotted the reaction rates from our calculations as a function of core 
density. In this log- log plot, rates are given as number of reactions per baryon per second, 
and the core density is given in terms of pio. As mentioned above, FCNC interactions are 
started in the simulation at a density of pio = 50. The dotted curve in the plot shows the 
unblocked electron capture reaction rate, i.e. the reaction rate computed as if there was 
no blocking of the final state z/g. This should not be confused with the actual net rate of 
electron capture (neutronization rate) in the core which is, of course, affected by blocking 
and by the reverse z/g capture reaction. We show the unblocked reaction rate to illustrate 
the relative size of the FCNC rates. Note that the FCNC rates for values of e < 10"'^ are 
smaller than the electron capture rate. Therefore, for values of e < 10~^, our assumption 
for computing AYe (where an electron capture is assumed to occur immediately when an 
electron neutrino changes flavor) is definitely justified. We expect that the assumption is 
also valid for values of 10^'^ < e < 10^^, though this range of epsilon is more complicated. 
For example, u^^j. neutrinos can undergo flavor changing scattering into u^s and sometimes 
fill holes before electron captures can. This scenario is discussed further in Section IIIII This 
range in e would be best modeled by a more detailed simulation which calculates neutrino 
transport. We anticipate that such a simulation would very our assumption. Finally, for 
values of e > 10^^, the FCNC cross section is comparable to the electron capture cross 
section and our assumption definitely does not hold. We therefore do not include these 
values of e in our simulation. 

From Figure [1] we see that even though the FCNCs are weaker and slower than Standard 
Model weak interactions, flavor changing scattering can nevertheless be significant. For 
e = 10~^, we see that at a density of pio = 100 there are ~ 10 flavor changing scatterings 
per baryon per second. At this density there are ~ 10^^ baryons per cubic centimeter in 
the core. Clearly then, there are a large number of electron neutrinos changing flavor as the 
core passes through this density. Of course, even greater numbers of neutrinos will change 
flavor as the collapse proceeds to higher densities. 
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III. FCNC-INDUCED REDUCTION IN ELECTRON FRACTION AND ALTER- 
ATIONS IN CORE PHYSICS 



Considering the FCNC rates in the core is useful in demonstrating that there are indeed 
a large number of neutrinos changing flavor, but the quantity which can be most important 
for the dynamics of the supernova model is the electron fraction. We have calculated the 
core's total AFg stemming from FCNCs as a function of the FCNC coupling constant e. As 
mentioned in the previous section, we find the total Al^g for each collapse simulation by 
summing Eq. ([8]) from each density step. The results are plotted in Fig. O Curves are 
shown for three different simulations, up to final densities of pio = 350, pio = 1140 and 
Pio = 3800. The range of the continuous parameter e is 10^^ to 10~^. The figure has vertical 
lines which show current values of experimental constraint on the epsilon parameter. The 
dotted line from e = 1 to e = 10^^ is included for ease in interpreting the figure, but we did 
not include these e values in our calculation. These values of epsilon are not strictly covered 
by the limiting case of e~ capture rates being faster than the z/g flavor changing rate. We 
discuss the case of large values of e below. 

In our simulation we assume a maximum trapped neutrino fraction of Y^, = 0.05 for each 
of the three flavors. Therefore, by Eq. ([7]), the maximum reduction possible for the electron 
fraction is AYe = —0.1. Of course, if one were to consider a different value for the maximum 
trapped neutrino fraction, the maximum AY^. would be different. As we will discuss below, 
and as marked on the figure, a dynamically significant reduction in Y^ can be as low as 
AYe = —0.02 [10|]. When the simulation runs to pio = 3800, our results show that the 
maximum reduction in Y^. occurs even for values of e as low as e = 10^'^. To put this in 
context, e = 10"'^ is coincident with the best current experimental constraint on the z/g <-> 
channel and is orders of magnitude smaller than the best current experimental constraint 
on z/g z/t-. For couplings e < 10""^, we see that reduction in Yf> due to FCNCs ceases to be 
significant. 

Figure [2] also has curves produced from simulations which were run to lower final densities. 
A reason for restricting our calculations to lower final densities is that our schematic liquid 
drop model equation of state is more reliable at lower density. However, since in these cases 
the FCNCs are not active for as long, there is less time for reduction in Ye to accumulate. 
The figure shows that even if FCNCs are active for only a short duration after trapping, 
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significant and/or maximal reduction of can occur for values of e allowed by current 
experimental bounds. 

In Tables Hi] and IIIII we show data from the simulation with active FCNCs. These tables 
can be compared to Table [I] which comes from a simulation without FCNC interactions. 
Table [ITl and Table [I] both go to the same final density of pio = 3800. As can be seen in 
Table [Tll Yuf, + reaches the maximum level (and maximum reduction in Y^ is obtained) 
at a density of pio = 2.87 x 10^. At this density the mean nucleus has a mass of A = 334. 
For the simulation without FCNCs included, the mean nucleus has a mass of A = 292 at 
this same density. The mean nucleus is larger in the simulation with FCNCs because the 
increased electron capture has caused nuclei to be more neutron rich. For the larger value 
of e = 10~^, Table [XTTl shows that maximum reduction in Ye occurs already at a density of 
pio ~ 200. At this density the mean nucleus has a less exotic size, A ^ 123. 

For our scenario where electron capture is fast compared to FCNCs, whenever an electron 
neutrino changes flavor, the hole in the z/g sea is assumed to be filled immediately by a 
produced via electron capture. We have ignored the possibility of mu or tau neutrinos 
undergoing FCNC scattering and changing into electron neutrinos, thus filling holes before 
electron capture can occur. It is most likely that holes would be filled by z/g's produced 
by electron capture. If a situation arises in the core where some holes are filled by mu or 
tau neutrinos that changed to electron neutrinos, the number of z/g's in the core would still 
remain the same. The z/g fraction still remains fixed at its maximum value of Y^^ = 0.05 and 
Eq. remains valid. In such a situation it is still possible for maximal reduction of Ye to 
occur. 

Situations other than those covered by our limiting case are also possible. These scenarios 
could occur for values of e close to 10~^. One example is if large FCNC rates cause significant 
numbers of mu and tau neutrinos to change flavor and seriously compete with electron 
capture in filling holes in the u^, sea. In this case, Y^^ still remains fixed at its maximum value 
and Eq. (JTj) is again valid, but maximal reduction in Yg would not occur. Another example 
is if the FCNC rates are greater than the electron capture rate. In this case we expect that 
the neutrino seas would equilibrate and all reach the same level, lowering the Fermi level of 
the z/g sea before electron capture reactions can replenish it. Instead of remaining fixed, Y^^ 
would be lowered initially. Maximum reduction in electron fraction may or may not occur. 
Equation ([7]) is not valid in this situation and the lepton number distribution in the core 
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would be different from that of our limiting case. Even if maximum possible reduction of 
does not occur during infall, we stress again that it only takes AYe ~ —0.02 to produce a 



significant alteration in core physics lOl 
According to the current modelj^, J s], 



mm 



for core collapse supernovae^eduction 



[il,Q,[3. 



of the core's electron fraction during infall will hinder the supernova explosionj 
Electrons in the core infiuence the collapse dynamics through the degeneracy pressure they 
provide. In particular, the number of electrons determines the size of the homologous core, 
Mhc ~ 5.8Y^Mq. a lower electron fraction, and consequentially smaller homologous core, 
hinders the explosion in two ways. If the inner core is smaller, there is more material in the 
outer region of the core for the shock to photodissociate before reaching the outer envelope 
of the star. Therefore, the shock has less energy available to eject the outer envelope and 
cause the explosion. A smaller inner core also has a smaller gravitational potential and so 
the outer core material has lower infall kinetic energy. The infall energy gets converted to 
the initial outgoing energy of the shock wave at bounce. (The initial shock energy scales 
as Y^^^^ as shown in Ref. 131].) Therefore, a smaller inner core results in a weaker shock. 
A weaker shock and greater loss of energy for the shock during its progression through the 
outer core may make an explosion more difficult to obtain. 

We have seen that including FCNCs in the supernova model causes Ye to be lowered 
and thus disfavors a successful explosion, or at least, can significantly alter the model. The 
FCNCs change the core's lepton number content. A standard collapse model would suggest 
that at bounce there would be a net electron lepton number in the core, but no net mu 
or tau lepton number. By contrast, with FCNCs there could arise significant net mu and 
tau lepton numbers resident in seas of z/^'s and z/^-'s. Of course, in this case we would still 
have sizeable electron lepton number residing in the electrons and the z/g's. This might have 
an interesting effect on the expected supernova neutrino signal, as speculated in Ref. j^. 



Neutrino medium-enhanced fiavor mixing 



oscillations) a 



30ve the neutron star will occur 



and will affect the signal [19|, l20|, UM, 122, l23|, \2S, l25|, l26|, \2m, l28| . Since neutrino mass-squared 
differences are known and mixing parameters (i.e., ^13) may be better constrained in the 
future, it may be possible to predict the effects of fiavor mixing and, upon detection of a 
supernova signal, subtract these out to identify signatures of FCNCs. An FCNC-engendered 
excess of z/^'s and/or z/^-'s might also result in altered neutrino energy/entropy transport in 
the proto-neutron star. 
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IV. DETAILED SUMMARY AND DISCUSSION OF APPROXIMATIONS 



We now summarize the approximations we have made and compare our calculation with 
modern detailed simulations. List of approximations: 

• One-zone 

• Constant collapse rate 

• Fermi Dirac distributions for neutrinos 

• Instantaneous trapping 

• Liquid drop model with mean representative nucleus 

• vacuum cross sections for neutrino nucleus interactions 

• ignoring medium effects and nuclear collective modes 

• electron capture rate is faster than flavor changing rate 

We first point out that even though our calculation has approximations, it is based on the 

n 

fact [1| that during the infall stage of collapse the core is a low entropy environment. Our 
equation of state is simple and based on Ref. The Ref. [l| EOS captures the essential 
features of more sophisticated treatments, for example Ref. 15|. Moreover, there is no 
difference in the underlying physical principals upon which our calculation and the detailed 
numerical simulations are based. 

Large detailed simulations follow the profile of the whole core, rather than a single zone 
as we have done. They calculate neutrino distributions and transport rather than imposing 
and modifying the distributions by hand. Neutrino trapping occurs "naturally" and the 
rate of the collapse is calculated self consistently within these simulations. In our case, we 
have not calculated these parameters in the presence of FCNCs. This is an approximation 
since we expect the physical effect of z/g u^ ^. will give some feedback on the system and 
these parameters will be modified. We point out however that the values we imposed (for 
example in the range 0.3-0.35, and collapse rate as in Eq. flA12l) are consistent with 
modern simulations. 
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As for nuclear matter, large detailed simulations treat this in the same way as we have. 
Nuclear matter at high density and temperature is still a wide area of research. Representing 
the nuclear matter in the core by an average nucleus, using a vacuum cross section for this 
nucleus which is in the core, and ignoring any other medium effects is an approximation. The 
large simulations make these same approximations when modeling the standard interactions. 
We have accounted for one medium effect: multiple target scattering when the neutrino's 
wavelength is comparable to the nuclear separation distance. 

To summarize, we have made approximations involving the thermodynamic profile and 
the core's dynamical quantities, neutrino transport, and nuclear matter. The thermody- 
namic and dynamical quantities we have used are based on the assumption of a low entropy 
core and are consistent with modern simulations. Our treatment of nuclear matter is similar 
to that in large scale numerical simulations. Neutrino transport near the trapping point 
is the issue that could be most improved upon by modeling FCNCs in a full simulation 
with Boltzmann neutrino transport. Finally, we point out that if a full simulation were 
performed to model these interactions, our qualitative conclusions would not change. Fur- 
thermore, our quantitative results are conservative and we expect a full simulation would 
indicate even smaller values of epsilon for which flavor changing significantly impacts core 
physics. This is because we have stopped our calculation at modest densities but larger 
numerical simulations would be able model these ultra high density regimes. 



V. CONCLUSIONS 



We have used a one- zone core collapse simulation to investigate some effects of including 
neutrino flavor changing interactions in the supernova model. We have calculated the re- 
duction in Ye as a function of the coupling constant e for collapse simulations that run up to 
density p = 3.8 x lO^^g/cm'^. For values of the interaction coupling constant e ^ 5 x 10~^ in 
either the z/e ^ or Uf. ^ Vj. channel we have found that maximal reduction in the core's 
electron fraction can occur. (See Fig. [2l) 

This work gives a more accurate and quantitative calculation of the effects of FCNCs 
than do the qualitative estimates given in Ref. P|. Here, we are able to account for the 
FCNC rate's dependence on density, and the feedback on the rates as the core becomes 
more neutron-rich as a result of increased net electron capture. However, a more accurate 
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treatment of these interactions is possible and is warranted. Some of our approximations 
were made for ease of calculation, while others were made to handle physics that is not 
yet well understood. Even with our conservative treatment (e.g., not following FCNCs 
beyond a density p = 3.8 x lO^^g/cm^), we were able to demonstrate how strong the effect 
of neutrino flavor changing interactions can be on the infall epoch physics. At very high 
densities near core bounce, we expect FCNCs still to be appreciable and to continue to cause 
net reduction in Yg. If a full simulation was preformed which included FCNCs and ran all 
the way to bounce density, properly accounting for neutrino scattering with nuclear matter, 
our results lead us to believe that significant reduction in Ye would occur for values of e even 
smaller than we have found here. 

The EOS and neutrino scattering cross sections in nuclear matter in the core are open 
areas of research. Some current simulations have more accurate treatments of these issues 
than we have used here. However, obtaining reliable cross sections for neutrino scattering 
with nuclear matter via Standard Model interactions remains problematic, in part because 
of the difficulty inherent in modeling nuclear matter. The standard neutrino interactions 
are treated with approximations, just as we have treated non-standard interactions with 
approximations. Our approximations are not a result of mysterious properties of FCNCs, 
but rather stem from uncertainties in matter at high density. The biggest uncertainty in 
our calculation does not come from the computational approximations in our model, but 
rather from this lack of knowledge. We point out these issues to differentiate physical 
approximations from computational approximations. When accurate and reliable EOS and 
compositions in nuclear matter in the core are available, standard and non-standard types 
of neutrino scattering can be correctly accounted for. 

A full supernova simulation, with neutrino transport and hydrodynamics, is needed to 
properly show all the effects of neutrino FCNCs. There are many pieces of known physics 
that are being tested for relevance in explaining supernova explosions 29|]. The supernova 
model cannot be used as a means of discovering or constraining new physics until known 
physics has been included and tested in simulations. Such simulations can treat neutrino 
trapping more realistically than we have done. By keeping track of neutrino distributions, 
such a simulation could handle the issue, discussed in Sec. Illlt of mu and tau neutrinos 
changing flavor and filling holes in the z/g sea before electron captures can occur. There 
is also a neutrino FCNC interaction with electrons j^. This additional opacity source for 
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neutrino flavor changing could be modeled easily in a full simulation. For all of these reasons, 
a better result for the reduction in Ye could be obtained. More sophisticated simulations 
also may reveal the fate of the shock, as well as changes to the thermodynamic proflle of the 
core. We used a constant collapse rate in our simulation, but in fact the pressure changes 
resulting from a continually decreasing would cause a non-uniform collapse rate. A full 
simulation would be able to follow the actual rate of collapse, and any consequences of a 
non-uniform collapse rate. Finally, a full simulation would provide neutrino specta which 
could reveal some signature of FCNCs in a supernova signal. The work presented in this 
paper will serve as a guide to preparing such a full simulation. 

Our results cannot be construed as either favoring or eliminating the existence of FCNCs. 
However, they do show that including FCNCs in the current supernova model could cause 
major changes to the model and its predictions. It is possible that data from a supernova 
signal could be used to constrain new physics such as FCNCs. On the other hand, new 
physics, such as what may be discovered at the LHC, might be required for successful 
explanation of supernovae. 
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APPENDIX A: ONE ZONE COLLAPSE PHYSICS 
1. Equation of State 

In the liquid drop model we can express the energy of a single nucleus as a sum of bulk, 

surface, and Coulomb terms, 

Wr,iYe, pn, Vn, u) = W^buik + W^surf A'/' + W^coui^'/'- (Al) 

Here, is the nuclear volume, u is the fraction of the total volume occupied by nuclei, A 
is the nuclear mass number, and VFsurf and VFcoui are coefficients of the surface and Coulomb 
energies, respectively. Deflning pjv as the density inside nuclei, we have A — pnVn and 
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u = p/pn- The coefficients W^urf and Wcoui are each functions of Ye, pn,Vn,u, and the 
number density of neutrons, We follow Ref.s Q, Q here and take these to be (in MeV) 



W^surf ~ 290 



i-x„ 



2 r 



Weoui ~ 0.75 (y^^) (1 - 0.236p}^' + 0.00194pi2), (A3) 

where pu = p/(10^^g/cm'^). The value for the mean nuclear mass is found from the free 
energy minimization condition, W^surf = 2VFcoui- The mean nuclear mass is 

A ^ 194 (y^)' (1 - 0.236p}/')-^ (A4) 

We follow Refs. Q, [3| and take the kinetic chemical potential (i.e., without rest mass) for 
neutrons to be (in MeV) 



Pn ~ -16 + 125 0.5 



1-X 



_ 125 (0.5-1 ^)) - CM) '-'^'l^-f^. (A5) 



X,J) V2/11/3/ 1 - (y;/(i - A'„)) 

In this expression we have neglected an additional term, — [W^surfM^^^(l — m^/^)]/[4A^/^(1 — 
3u^/^/2 + u/2)]. This is justified when u is small. Neglect of this term will cause some 
inaccuracy at the highest densities shown in our tables. The neutron-proton kinetic chemical 
potential difference is (in MeV) 



p = 250 0.5 



l-Xr, 



Ye V'^f Ye \-M-2(y;/(i-x„)) 



+ 2 



,i-x„; vi-x„y i-(n/(i-x„)) 

The mass fraction of free neutrons in the dilute limit is 



(A6) 



T3/2 

Xn ^ 79 e^"/^, (A7) 

Pio 

where pio = p/(10^°g/cm^). At very high densities neutron degeneracy becomes important 
and this expression will be inadequate. Likewise, in the dilute limit the free proton mass 
fraction is 

Xp ^ X^e^/^. (A8) 
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2. Collapse Rate 



.0 



Following Refs. [l| and [iSj] we choose a collapse rate which is a fraction of the free fall 
rate. The free fall rate for a core with mass M interior to a radius R is 



R 
R 



where G is Newton's Constant. With p = M/{AttR^/3) we have 

dlnR 1 dlnp 



dt 



3 rft 



and so 



^.(224s-)(p,o)^/^ 



(A9) 



(AlO) 



(All) 



However, as discussed in Ref. ISj, the actual collapse rate is smaller than this. This is 
because degenerate electrons provide pressure which slows the collapse. Hence, we again 
follow Ref. [13|] and take 

:iOOs-i)(pio)'/' (A12) 



dlnp 



dt 

for the collapse rate. This collapse rate was chosen in Ref. m to represent the rate found 

h n n 

from then existing numerical simulations. Contemporary simulaions [7|, l8|, l9[ give similar 
collapse rates. The gross features of collapse are set by the low entropy conditions and 
the rate of electron capture on heavy nuclei (a quantity dominated by the high electron 
degeneracy) . 



3. Entropy 



The change in entropy per baryon with density, derived in ISj, is adopted here: 



(dS\ (-dYAI . (l-.Y„-AgAB 

+ 1,. _ ^ _ _ iX^IX.)] ^^^^^^^^^^:'l_^^^^^ \ (A13) 

The neutron-proton mass difference is (5m„p 1.293 MeV, and Afp and AJ^ are neutrino 
energy loss rates for electron capture on free protons and heavy nuclei, respectively. The 
total (with rest mass) electron chemical potential pe is given by 

Pe~ll.lMeV(pioye)'/'. (AM) 
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The entropy per baryon (in units of Boltzmann's constant k) for the system of nuclei, 
nucleons, and electrons employed in our one- zone collapse is 



S 



+ 0.2467(1 - Xn 



T 



+ 



+ 



2.5(1 (1-XO, 
- H : in 



A 



2.5X„ + X„ln 



A 

79.07(T)3/^ 



39.49^5/2 (T)3/2 

Plo(l - Xn) 



(A15) 



where T and /ig are in MeV. 



APPENDIX B: NEUTRINO NUCLEUS SCATTERING IN MEDIUM 

The general FCNC Lagrangian is given in Eq. ([T]). For illustrative purposes, we consider 
the u — d term from this Lagrangian and calculate the FCNC cross section for this channel. 
For the sake of generality, we drop subscripts on the parameter e. The differential cross 
section for flavor changing neutrino-nucleus coherent scattering is then 

-''^^{2N + ZfEl{l + cose)F\q). (Bl) 



dcOs9 TT 

Here, E,^ is the incident neutrino energy, and Z are the numbers of neutrons and protons, 
respectively, in the nucleus, 6 is the scattering angle, and q is the momentum transfer given 
by 

q = V2E^{l-cosey^^. (B2) 

In Eq. (IBip . F{q) is a form factor. In the calculations done in Ref. j^, this form factor was 
set to unity. 

In medium, when the neutrino wavelength is comparable to the separation distance be- 
tween nuclei, that is l/E^, ^ R^, interference from multiple nucleus scattering can be ac- 
counted for by modifying the form factor. This is done by subtracting a term from the 
original form factor [l8|, 

F{q) = F{q) - , M<lRc) - {^Rc)cosiqR.) 

and then replacing F'^{q) by F'^{q) in the differential cross section. In our case, with F{q) ^ 
1, we have 

,sin(g_Rc) — {qRc)cos{qR, 



Jl^ = f^^2N + Z)'Elil + cose)\l-3-' ' 
dcosfc' vr 



{qRcY 



(B4) 
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The cross section is then given by 

a=t^{2N + ZfElX, (B5) 

TT 

where 1 is the integral of the kinematical factor (1 +cos^) multiphed by F'^{q). Using Eq. 
(IB2I) to change variables, the integral is 

The integral X is evaluated numerically at each density step using that step's values for Ey 
and Rc = -Dsep/2. 

To illustrate how the correction for multiple target scattering can change the cross section 
we have included plots of the correction for particular values of core density. In Figure [3] we 
have plotted the corrected form factor when the core is at a density pio = 748. From Table 
[Iwe see that = 35.10 MeV and Rc = 20.56 fm = 0.112 MeV~^ at this particular density. 
From Figure [3] we clearly see that for small scattering angles (low values of q) the interference 
is destructive and the differential cross section is reduced. However, we must integrate over 
q. In Figure m we have plotted the kinematical factor, {2q/ E^ — q^ /2E^) multiplied first by 
the uncorrected form factor F'^{q) = 1 and then by the corrected form factor F'^{q). This 
plot was made for values of E^, and Rc corresponding to density pio = 125. Table [T] gives 
Ey = 17.83 MeV and Rc = 33.01 fm = 0.168 MeV"^ at pio = 125. At this lower density we 
see that the differential cross is lower over the whole range of q. The interference causes the 
integrated cross section to be reduced. 
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H10 

FIG. 1: (color online). Neutrino FCNC scattering rates as a function of density for the indicated 
values of FCNC coupling e. The solid curve shows the unblocked electron capture reaction rate. 
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FIG. 2: (color online). Magnitude of the net decrease in electron fraction (— AYJ,) as a function of 
e, for collapse up to various final densities. Curves become dotted for e > 10^^ because we have 
not included this range of e in our calculation; the dotted line is included for ease of reading the 
figure. The horizontal line indicates a threshold change in electron fraction {—AYe > 0.02) beyond 
which significant alteration in core physics can be expected. The two vertical lines give the current 
experimental bounds on e = eer (for Ve i^r) and e = (for z/e f^). 
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FIG. 3: (color online). Corrected form factor as a function of momentum transfer q for core density 
pio = 748. At this density Rc = 20.56 fm, and = 35.10 MeV. 
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FIG. 4: (color online). The solid curve shows kinematical factor multiplied by uncorrected form 
factor and the dotted curve shows kinematical factor multiplied by the corrected form factor, 
each plotted as functions of momentum transfer g, for core density piQ = 125. At this density 
Rc = 33.01 fm, and = 17.83 MeV. 
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PlO 


Ye 


Yu.^ 




T (MeV) 


s/k 




/Un (MeV) 


A 


Tnuc (fm) 


L>sep (fm) 


(McV) 


K (fm) 


0.37 


0.4200 


0.0000 


0.0000 


0.660 


0.90 


0.0024 


-7.10 


67. 


4.07 


387.26 


0.00 


inf 


0.58 


0.4200 


0.0000 


0.0000 


0.720 


0.90 


0.0041 


-7.16 


68. 


4.08 


334.17 


0.00 


inf 


0.91 


0.4200 


0.0000 


0.0000 


0.782 


0.90 


0.0059 


-7.22 


68. 


4.08 


288.44 


0.00 


inf 


1.42 


0.4200 


0.0000 


0.0000 


0.847 


0.91 


0.0080 


-7.29 


68. 


4.09 


249.08 


0.00 


inf 


2.22 


0.4198 


0.0000 


0.0000 


0.914 


0.91 


0.0101 


-7.35 


69. 


4.10 


215.22 


0.00 


inf 


3.47 


0.4193 


0.0000 


0.0000 


0.984 


0.91 


0.0123 


-7.38 


70. 


4.11 


186.15 


0.00 


inf 


5.43 


0.4177 


0.0000 


0.0000 


1.055 


0.91 


0.0149 


-7.34 


71. 


4.13 


161.29 


0.00 


inf 


8.50 


0.4143 


0.0000 


0.0000 


1.124 


0.91 


0.0185 


-7.19 


72. 


4.17 


140.15 


0.00 


inf 


13.30 


0.4084 


0.0000 


0.0000 


1.190 


0.92 


0.0238 


-6.88 


75. 


4.21 


122.23 


0.00 


inf 


20.81 


0.3998 


0.0000 


0.0000 


1.251 


0.94 


0.0313 


-6.42 


78. 


4.27 


107.02 


0.00 


inf 


32.57 


0.3884 


0.0000 


0.0000 


1.313 


0.97 


0.0425 


-5.84 


82. 


4.34 


94.10 


0.00 


inf 


50.96 


0.3736 


0.0009 


0.0000 


1.375 


1.01 


0.0591 


-5.15 


87. 


4.43 


83.18 


3.73 


331.39 


79.74 


0.3559 


0.0185 


0.0000 


1.372 


1.01 


0.0757 


-4.18 


94. 


4.54 


73.96 


11.95 


103.52 


124.77 


0.3351 


0.0394 


0.0000 


1.351 


1.01 


0.0990 


-3.12 


103. 


4.68 


66.17 


17.83 


69.37 


195.24 


0.3244 


0.0501 


0.0000 


1.430 


1.01 


0.1130 


-2.59 


111. 


4.80 


58.78 


22.43 


55.15 


305.51 


0.3244 


0.0501 


0.0000 


1.622 


1.01 


0.1120 


-2.53 


119. 


4.92 


51.82 


26.04 


47.50 


478.07 


0.3244 


0.0501 


0.0000 


1.834 


1.01 


0.1096 


-2.42 


130. 


5.07 


45.96 


30.23 


40.91 


748.07 


0.3244 


0.0501 


0.0000 


2.064 


1.01 


0.1054 


-2.25 


146. 


5.27 


41.11 


35.10 


35.24 


1170.57 


0.3244 


0.0501 


0.0000 


2.310 


1.01 


0.0995 


-2.00 


171. 


5.55 


37.22 


40.75 


30.36 


1831.70 


0.3244 


0.0501 


0.0000 


2.564 


1.01 


0.0917 


-1.69 


212. 


5.96 


34.34 


47.31 


26.15 


2866.21 


0.3244 


0.0501 


0.0000 


2.821 


1.01 


0.0825 


-1.30 


292. 


6.63 


32.79 


54.92 


22.52 



TABLE I: Sample calculation with no flavor changing coupling, i.e., e = 0. Values are density pio, 
electron fraction Y^, mu plus tau neutrino fraction Y^^ + Y^,^, temperature T in MeV, entropy per 
baryon s/k \n units of Boltzmann's constant, free neutron fraction X„, neutron kinetic chemical 
potential /7,„ in McV, mean nuclear mass A, mean nuclear radius rnuc m fm, separation distance 
between nuclei -Dsep m fm, average neutrino energy E^, in MeV, and average neutrino de Broglie 
wavelength Xi, in fm. 
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PlO 


Ye 






T (MeV) 


s/k 




/i„ (MeV) 


A 


Tnuc (fm) 


Dsep (fm) 


(MeV) 


K (fm) 


0.37 


0.4200 


0.0000 


0.0000 


0.660 


0.90 


0.0024 


-7.10 


67. 


4.07 


387.26 


0.00 


inf 


0.58 


0.4200 


0.0000 


0.0000 


0.720 


0.90 


0.0041 


-7.16 


68. 


4.08 


334.17 


0.00 


inf 


0.91 


0.4200 


0.0000 


0.0000 


0.782 


0.90 


0.0059 


-7.22 


68. 


4.08 


288.44 


0.00 


inf 


1.42 


0.4200 


0.0000 


0.0000 


0.847 


0.91 


0.0080 


-7.29 


68. 


4.09 


249.08 


0.00 


inf 


2.22 


0.4198 


0.0000 


0.0000 


0.914 


0.91 


0.0101 


-7.35 


69. 


4.10 


215.22 


0.00 


inf 


3.47 


0.4193 


0.0000 


0.0000 


0.984 


0.91 


0.0123 


-7.38 


70. 


4.11 


186.15 


0.00 


inf 


5.43 


0.4177 


0.0000 


0.0000 


1.055 


0.91 


0.0149 


-7.34 


71. 


4.13 


161.29 


0.00 


inf 


8.50 


0.4143 


0.0000 


0.0000 


1.124 


0.91 


0.0185 


-7.19 


72. 


4.17 


140.15 


0.00 


inf 


13.30 


0.4084 


0.0000 


0.0000 


1.190 


0.92 


0.0238 


-6.88 


75. 


4.21 


122.23 


0.00 


inf 


20.81 


0.3998 


0.0000 


0.0000 


1.251 


0.94 


0.0313 


-6.42 


78. 


4.27 


107.02 


0.00 


inf 


32.57 


0.3884 


0.0000 


0.0000 


1.313 


0.97 


0.0425 


-5.84 


82. 


4.34 


94.10 


0.00 


inf 


50.96 


0.3736 


0.0009 


0.0000 


1.375 


1.01 


0.0591 


-5.15 


87. 


4.43 


83.18 


3.73 


331.39 


79.74 


0.3559 


0.0185 


0.0000 


1.372 


1.01 


0.0757 


-4.18 


94. 


4.54 


73.96 


11.95 


103.52 


124.77 


0.3349 


0.0393 


0.0002 


1.349 


1.01 


0.0992 


-3.11 


103. 


4.68 


66.18 


17.82 


69.39 


195.24 


0.3230 


0.0501 


0.0013 


1.419 


1.01 


0.1148 


-2.53 


111. 


4.81 


58.86 


22.43 


55.14 


305.51 


0.3207 


0.0501 


0.0036 


1.593 


1.01 


0.1170 


-2.38 


120. 


4.93 


52.03 


26.04 


47.49 


478.07 


0.3165 


0.0501 


0.0079 


1.774 


1.01 


0.1208 


-2.08 


132. 


5.09 


46.37 


30.24 


40.91 


748.07 


0.3080 


0.0501 


0.0163 


1.951 


1.01 


0.1296 


-1.56 


150. 


5.32 


41.86 


35.10 


35.24 


1170.57 


0.2906 


0.0501 


0.0338 


2.110 


1.01 


0.1528 


-0.64 


180. 


5.65 


38.66 


40.75 


30.35 


1831.70 


0.2526 


0.0501 


0.0717 


2.269 


1.01 


0.2201 


0.91 


235. 


6.17 


37.36 


47.31 


26.14 


2866.21 


0.2242 


0.0501 


0.1001 


2.612 


1.01 


0.2724 


2.22 


334. 


6.94 


37.06 


54.93 


22.52 



TABLE II: Same as Table H] but now with flavor changing coupling e = 10 
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PlO 


Ye 


Yu. 




T (MeV) 


s/k 




(MeV) 


A 


Tjiuc (fm) 


L>scp (fm) 


(MeV) 


K (fm) 


0.37 


0.4200 


0.0000 


0.0000 


0.660 


0.90 


0.0024 


-7.10 


67. 


4.07 


387.26 


0.00 


inf 


0.58 


0.4200 


0.0000 


0.0000 


0.720 


0.90 


0.0041 


-7.16 


68. 


4.08 


334.17 


0.00 


inf 


0.91 


0.4200 


0.0000 


0.0000 


0.782 


0.90 


0.0059 


-7.22 


68. 


4.08 


288.44 


0.00 


inf 


1.42 


0.4200 


0.0000 


0.0000 


0.847 


0.91 


0.0080 


-7.29 


68. 


4.09 


249.08 


0.00 


inf 


2.22 


0.4198 


0.0000 


0.0000 


0.914 


0.91 


0.0101 


-7.35 


69. 


4.10 


215.22 


0.00 


inf 


3.47 


0.4193 


0.0000 


0.0000 


0.984 


0.91 


0.0123 


-7.38 


70. 


4.11 


186.15 


0.00 


inf 


5.43 


0.4177 


0.0000 


0.0000 


1.055 


0.91 


0.0149 


-7.34 


71. 


4.13 


161.29 


0.00 


inf 


8.50 


0.4143 


0.0000 


0.0000 


1.124 


0.91 


0.0185 


-7.19 


72. 


4.17 


140.15 


0.00 


inf 


13.30 


0.4084 


0.0000 


0.0000 


1.190 


0.92 


0.0238 


-6.88 


75. 


4.21 


122.23 


0.00 


inf 


20.81 


0.3998 


0.0000 


0.0000 


1.251 


0.94 


0.0313 


-6.42 


78. 


4.27 


107.02 


0.00 


inf 


32.57 


0.3884 


0.0000 


0.0000 


1.313 


0.97 


0.0425 


-5.84 


82. 


4.34 


94.10 


0.00 


inf 


50.96 


0.3736 


0.0009 


0.0000 


1.375 


1.01 


0.0591 


-5.15 


87. 


4.43 


83.18 


3.73 


331.39 


79.74 


0.3550 


0.0184 


0.0010 


1.365 


1.01 


0.0766 


-4.13 


94. 


4.55 


74.03 


11.92 


103.77 


124.77 


0.3192 


0.0355 


0.0197 


1.215 


1.01 


0.1197 


-2.38 


106. 


4.73 


67.34 


17.23 


71.78 


195.24 


0.2313 


0.0429 


0.1003 


0.811 


1.01 


0.3032 


0.02 


123. 


4.97 


65.94 


21.30 


58.08 


305.51 


0.2265 


0.0477 


0.1003 


1.001 


1.01 


0.3138 


0.19 


132. 


5.10 


58.53 


25.61 


48.29 


478.07 


0.2241 


0.0500 


0.1003 


1.239 


1.01 


0.3165 


0.41 


145. 


5.26 


52.10 


30.22 


40.93 


748.07 


0.2241 


0.0500 


0.1003 


1.528 


1.01 


0.3104 


0.68 


164. 


5.48 


46.58 


35.09 


35.25 



TABLE III: Same as Table H] but now with flavor changing coupling e = 10 
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